Cavity Quantum Electrodynamics with a Rydberg blocked atomic ensemble 



Christine Guerlin 1 ' 4 , Etienne Brion 2 , Tilman Esslinger 1 and Klaus M0lmer 3 
1 Institute for Quantum Electronics, ETH Zurich, CH-8093 Zurich, Switzerland, 
2 Laboratoire Aime Cotton, CNRS, Campus a 1 ' Or say, 91405, Or say, France, 
3 Lundbeck Foundation Theoretical Center for Quantum System Research, 
Department of Physics and Astronomy, University of Aarhus, DK-8000, Denmark, 
A Thales Research and Technology, Campus Poly technique, 
1 avenue Augustin Fresnel, 91767 Palaiseau Cedex, France. 
(Dated: June 21, 2010) 

We propose to implement the Jaynes-Cummings model by coupling a few-micrometer large atomic 
ensemble to a quantized cavity mode and classical laser fields. A two-photon transition resonantly 
couples the single-atom ground state \g) to a Rydberg state \e) via a non-resonant intermediate state 
\i), but due to the interaction between Rydberg atoms only a single atom can be resonantly excited 
in the ensemble. This restricts the state space of the ensemble to the collective ground state \G) and 
the collectively excited state \E) with a single Rydberg excitation distributed evenly on all atoms. 
The collectively enhanced coupling of all atoms to the cavity field with coherent coupling strengths 
which are much larger than the decay rates in the system leads to the strong coupling regime of the 
resulting effective Jaynes-Cummings model. We use numerical simulations to show that the cavity 
transmission can be used to reveal detailed properties of the Jaynes-Cummings ladder of excited 
states, and that the atomic nonlinearity gives rise to highly non-trivial photon emission from the 
cavity. Finally, we suggest that the absence of interactions between remote Rydberg atoms may, 
due to a combinatorial effect, induce a cavity- assisted excitation blockade whose range is larger than 
the typical Rydberg dipole-dipole interaction length. 

PACS numbers: 32.80.Ee, 32.80. Qk, 37.30.+i, 42.50.Pq 



I. INTRODUCTION 

The Jaynes-Cummings model (JCM) pQ provides the 
general framework to describe the interaction of a two- 
level system - such as an atom, with a quantum har- 
monic oscillator - e.g. a quantized cavity mode in cavity 
quantum electrodynamics (CQED). It allows to explain 
very specific behaviors of such systems as, for example, 
the collapses and revivals of Rabi oscillations [2j [3] . The 
JCM was also used to describe various physical situations 
outside the field of CQED, such as the coupling of the 
internal states to the centre-of-mass vibrational levels of 
a trapped ion subject to a laser beam in the Lamb-Dicke 
regime (for a review, see [4 ), or, more recently, a laser- 
driven electron floating on liquid Helium [5 . Precisely 
because of this universality, the JCM also appears now 
as one of the key ingredients for applications in quantum 
information processing. 

While many successful CQED experiments have been 
carried out with microwave cavities, resonant with tran- 
sitions between atomic Rydberg excited states [6[ [7] and 
with superconducting "artificial atoms" [SJE], it is harder 
to reach the strong coupling regime with single atoms and 
optical cavities. This regime is, however, of considerable 
practical interest, as it makes it possible for the cavity 
to work as an interface between the optical photons, fly- 
ing qubits, and the atomic states, stationary qubits, in 
quantum computing and communication networks [ID] . 

A viable path to obtain a highly coherent system in the 
optical regime is to use a collection of N atoms whose 
coupling strength to an optical field mode is magnified 
by the factor relative to the single-atom case [TT1 - 



[21]. For large TV, the collective excitation degree of free- 
dom of an atomic ensemble is approximately equivalent 
to a quantum harmonic oscillator and the system atoms- 
cavity therefore tends to be well-described by a simple 
quadratic Hamiltonian in the raising and lowering op- 
erators for the atomic and field excitations. Concepts 
to squeeze and entangle the atomic and field degrees of 
freedom by this Hamiltonian have been developed but 
the systems stay within the restricted family of Gaussian 
states [22] 123]. 

In this paper, we consider an atomic ensemble placed 
in an optical high-finesse cavity and investigate the non- 
linearity introduced by the interactions between atoms 
resonantly coupled to a high-lying Rydberg state. In 
relatively small ensembles, dipole-dipole interactions sig- 
nificantly shift the energy of states with two or more 
Rydberg excitations [23] . When states with two excited 
atoms are shifted far away from the optical resonance 
we observe the so-called Rydberg Blockade phenomenon 
[25H3T] , where the atomic ensemble effectively behaves as 
a single two-level system. The entire atoms-cavity sys- 
tem is therefore well-described by the Jaynes-Cummings 
model, and the collective enhancement factor com- 
pared to the single atom case may then lead to the strong 
coupling regime. 

The analyses we shall provide in this manuscript re- 
fer to a specific experimental implementation with real- 
istic physical parameters specified in Sec. II. The ex- 
pected signature of the strongly coupled ensemble on 
the optical transmission of the cavity is discussed and 
investigated numerically by Monte- Carlo simulations in 
Sec. III. In particular, we show that both the average 
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signal and the fluctuations present interesting, and per- 
haps surprising, features linked with the eigenstates of 
the Jaynes-Cummings model. In Sec. IV, we suggest 
a new mechanism by which the cavity can extend the 
range of the Rydberg blockade beyond the dipole-dipole 
interaction length. In Sec. V, we conclude and discuss 
briefly the vast range of possible investigations that can 
be made with the system and a few possible extensions 
of the theoretical model. 



II. PHYSICAL IMPLEMENTATION 

The principle of our proposed experiment is depicted 
in Fig. [I] A Bose-Einstein condensate consisting of 
N > 1000 87 Rb atoms is placed in an ultra-high finesse 
cavity [17] and is transversely illuminated by a homoge- 
neous classical laser field. The 5Si/ 2 ground atomic state 
denoted by \g) is assumed to be resonantly coupled to a 
Rydberg level \e) through a two-photon process, via the 
5P3/2 intermediate state \i). The transition \g) — >> of 
frequency uj 9 ^h, is non-resonantly coupled with the cou- 
pling strength go to a single quantized cavity mode with 
annihilation operator a and frequency uj c . The cavity 
mode is detuned by the amount A = uj c —uj g ^i while the 
transition \i) — >> |e) of frequency cj^ e , is non-resonantly 
driven by the laser field with the Rabi frequency ft and 
frequency uj\ detuned by the amount uj\ — 0Ji^ e = —A. 
Omitting at first the interatomic interactions, one can 
describe the physical situation by the following Hamilto- 
nian 
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written in the interaction picture with respect to Ho = 

hhj c {a)a + Y^=i + K^c+ui) Ylf=i \ e j)( e jl and 

in the rotating wave approximation. Dispersive shifts 
on both transitions have been neglected for simplicity. 
Assuming a large detuning from the intermediate state, 
A ^> T, where T _1 is the lifetime of we can adia- 
batically eliminate the unpopulated intermediate state, 
which leads to the effective Hamiltonian 



i^eff = h^a^\e j ){g j \+h.c. 
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where the collective mode atomic excitation is described 

by5t = ^Ef =1 |e J )(5,|. 

Eq. (1) is written under the assumption that all atoms 
have the same coupling coefficients go, Q to the laser 



and cavity field. This describes the situation of a Bose- 
Einstein condensate, when neglecting the photon recoil. 
Alternatively, one can consider a real space where the 
coupling strengths depends on the value of the mode 
functions at the location of the individual atoms, and g^ 
should be replaced by different coefficients g*£lj on each 
term in the sum in (2). It turns out, however, that in the 
limit of many atoms, one can define a collective operator 
like which is a weighted sum of the individual atomic 
raising operators, and which also obeys the oscillator like 
commutator relations to a good approximation. As long 
as the atoms do not move appreciably and change their 
coupling strengths on the time scale of interest for the 
experiment, this weighted collective atomic mode plays 
the same role as the ideal symmetric mode and its cou- 
pling is enhanced by the same factor y/~N and an extra 
mode dependent factor of order unity, see, e.g., [33] . 
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For low Rydberg excitation numbers, 

the Hamiltonian Eq. ([2| approximately describes the 
beam splitter coupling of two degenerate oscillators. 

The dipole-dipole interactions we have omitted so far 
are strong only between neighbouring Rydberg atoms. 
Their effect on an atomic sample whose size is at most of 
the order of the typical Rydberg-Rydberg dipole-dipole 
interaction length £r is to considerably shift the energy of 
multiply Rydberg excited states. Driving the transition 
\g) — >• \e) as described above, one can thus resonantly 
couple the collective ground state \G) = \gi, . . . ,gw) to 
the collective symmetric state with a single Rydberg ex- 
citation \E) = \G), only, as higher excited levels are 
too far detuned. This constitutes the Rydberg blockade 
phenomenon [25j [26]. Restricting the Hamiltonian Eq. 
([2| to the physically relevant atomic subspace {\G) , \E)} 
we get H'eff = %eff^ \E) (G\ + h.c. which corresponds 
to the interaction part of the Jaynes-Cummings Hamil- 
tonian describing the resonant interaction of the cavity 
field with a fictitious two- level "super atom" [32] . with 
the effective coupling constant g e R = V^gg^/A. We 
finally pass to a new rotating frame defined by the uni- 
tary state transformation exp {-itu c (a f a+ \E) (E\)} to 
obtain the full Jaynes-Cummings Hamiltonian 



Hjc = hhj c a [ a J r huj c \E){E\ 
+ [% eff a|^) <G|+ft.c] 



(3) 



with an effectively equal excitation energy of the field 
and ensemble atomic states. 

The parameters of this Jaynes-Cummings model imple- 
mentation can be tuned over a wide range. The collective 
coupling on the first transition \fN go scales indeed as the 
square root of the atom number and can thus be varied 
in the experiment up to the GHz range [17] . One possible 
scheme would be to excite the 70s Rydberg state from a 
1000 atoms BEC, with a cavity-intermediate state detun- 
ing of A = 2ir x 900 MHz, a single atom resonant coupling 
on that transition of go = 2ir x 10 MHz and a blue laser 
Rabi coupling Q = 2tt x 30 MHz. Excitation of the inter- 
mediate atomic state \i) can be neglected here since the 
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cavity detuning A is large compared to both the atomic 
radiative decay rate r = 2tt x 3 MHz from the intermedi- 
ate state and the one-photon collective coupling strength 
y/~N go ~ 2tt x 320 MHz. The resulting two-photon cou- 
pling g e ff ~ 2ir x 10 MHz is higher than both decay rates 
of the cavity k = 2tt x 1.3 MHz and of the Rydberg 
state 7 = 2tt x 0.55 kHz. This system reaches the strong 
coupling regime on the ground-to-Rydberg state two- 
photon transition. Within the BEC Thomas-Fermi ra- 
dius Vtf ~ 2 fim corresponding to an isotropic trapping 
frequency 2tt x 180 Hz, the minimum Rydberg interaction 
shift over the atomic sample is Ari ~ 2tt x 200 MHz, 
which spectrally separates the states with two or more 
Rydberg excited atoms from the singly excited ones. 



III. THE DRIVEN JAYNES-CUMMINGS 
MODEL 

In the previous section, we arrived at an effective de- 
scription of the ensemble-light interaction in terms of a 
two-level super atom whose Hilbert space consists of the 
collective ground and singly Rydberg excited states of 
the ensemble |G) , \E). The associated effective Jaynes- 
Cummings Hamiltonian Eq. Q is block diagonal, cou- 
pling only pairs of states (|G, n) , \E, n — 1)) with the 
same total number n of either photonic or atomic exci- 
tations. The frequencies uj n ,± of the dressed eigenstates 
|n, ±), obtained by direct diagonalization of Hjc, are 
given by 

u Ui ± = nuj c ± g e fiVn. (4) 

The resulting nonlinearity is depicted on Fig. [2] 

This system can be driven by a classical probe field, 
which either illuminates the atoms on the g — i transition 
or feeds the cavity via one of the cavity mirrors. In the 
following we shall consider the latter solution. Due to 
the cavity and atomic decay, the combined field-atom 
system will reach a steady state after a few relaxation 
times, i.e. a few microseconds for the specific system 
under consideration. If the driving field is weak, this state 
only slightly differs from the ground state |G, n = 0) and 
depending on the driving frequency, we may estimate the 
probability amplitude with which the system populates 
the excited eigenstates \n, ±) of the unperturbed Jaynes- 
Cummings model. 



A. One-photon resonance 

To excite the system, the driving frequency should 
match one of the two dressed state components at 
uj c ± g e ff? shifted by the so-called vacuum Rabi split- 
ting. Directly measuring the absorption or transmission 
of the cavity as a function of the driving frequency , 
thus reveals the collectively enhanced coupling strength 



g e ff. The width of the resonances is mainly governed 
by the cavity decay (as noted above the atomic decay 
rate 7 is negligible compared to the cavity decay rate 
k). With # e ff ~ 2tt x 10 MHz and k ~ 2n x 1 MHz, 
the two resonance lines are then clearly split, and it is 
possible to selectively excite one of the dressed states 
|±, n = 1) = (|G, 1) ± \E, 0)) /V2. If the probe field cou- 
ples to the cavity mode according to the Hamiltonian 

V r = a*e <wt a + ae- <w *a t , (5) 

the coupling strength of the ground state to either 
of the dressed states is the same, given by j3\ = 
(±,n = 1| V |G, 0) = a*/y/2. We therefore expect these 
states to be excited with a probability pi ~ \fii\ 2 /(S 2 + 
ft 2 /4), where S denotes the probe frequency detuning 
with respect to the dressed state eigenfrequency. The 
dressed state populated is a superposition of the atomic 
and field excited state, and the mean photon number in- 
side the cavity is expected to be (n) = \pi- 

To analyze the problem theoretically, we have carried 
out numerical Monte Carlo simulations of the dynam- 
ics of the driven atom+cavity system. The cavity de- 
cay occurs by emission of photons [34], and we simu- 
late the detection of these emission events by a photon 
counter in the Monte Carlo Wave Function (MCWF) for- 
malism [35, 36 j. Such simulations do not converge to a 
steady state, they rather present a stochastic dynamics 
with click events followed by transient evolution until the 
next click event is detected. When averaged over many 
independent realizations of this stochastic process, one 
recovers the predictions of the master equation [34] [37] . 
Moreover, the individual stochastic simulations are also 
useful, as they provide typical records of the randomly 
selected detection events which is, indeed, information 
similar to the one obtained in a real transmission exper- 
iment. 

Fig. [3] shows the result of a simulation, where we plot 
the mean photon number in the cavity (n), as a func- 
tion of time. The graph shows characteristic oscillations 
at the Rabi frequency of the transition, interrupted by 
sudden jumps corresponding to the detection of a pho- 
ton. At jumps the intra-cavity photon number drops to 
nearly zero, while the state of the system, in the dressed 
state superposition of |G, 1) and \E : 0) just before the 
measurement, collapses to the ground state |G, 0). One 
also observes an immediate consequence of these jumps: 
the system is unable to emit another photon until it has 
been reexcited, hence the clicks on the detector are pre- 
dicted to arrive at a rate given by the mean excited state 
population and to be antibunched on the time scale of 

B. Two-photon resonance 

Recently [38], Schuster et al. studied the Jaynes- 
Cummings model of a single atom in an optical cavity, 
and they identified a peak in the transmitted power, 
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FIG. 1: Left: A Bose-Einstein condensate inside an ultra- high finesse optical cavity is coupled in a two photon scheme to a 
Rydberg state, with the first photon being provided by a single quantized cavity mode and the second photon by a classical 
coupling beam. The system is interrogated by monitoring the cavity transmission of a weak probe beam on a photon counter. 
Middle: The quantized cavity mode has a frequency z/i and couples the collective atomic ground state \G) to an intermediate 
level | J) with a detuning A. The collective vacuum Rabi coupling is given by go\ /r N. A blue photon of frequency v<i drives 
the excitation to a Rydberg state \E) with a classical Rabi coupling Q. Right: The interaction between two Rydberg atoms 
leads to an energy shifted and broadened state \E) containing two excitations. This allows selectively addressing the single 
excitation state. The resulting situation can be described by a Jaynes-Cummings Hamiltonian. 



when the cavity is driven at one of the frequencies 
u c ± g e f? I \/2. Indeed, the driving field is then resonant - 
by two-photon absorption - with one of the dressed states 
\n = 2, zb) with respective energies UJ2,± = 2co> c ± g e ff V^2- 
This is a difficult experiment, both because this reso- 
nance is only a few decay widths away from the single 
photon resonance, and because the light atom interaction 
may cause heating of the atomic motion in the frequency 
range giving the optimal resolution of the resonance. 

With our parameters, the splitting is larger, and 
due to the long lifetime of the atomic Rydberg ex- 
cited states, the light induced heating will be much re- 
duced. The second-order transition between the ground 
state |G,0) and the dressed states |±,n = 2) is detuned 
from the intermediate states |±,n = 1) by the amount 
u c ± g ei ?V2 - (u c ±g ei ?/V2) = ±g eS (l - l/y/2). The 
two-photon coupling strength is hence estimated to be 
ft ~ ±(a/V2)a(l + >/2)/(0eff (1 - 1/v^)) ~ ±3a 2 /# eff . 
Accordingly, the steady state population of the doubly 
excited dressed state is p2 ~ |02| 2 /(^ 2 + ^ /2 ), where k' 
denotes the coherence decay rate of the photonic com- 
ponents of the dressed states. Since the dressed state is 
composed of equal weight components with a single and 
two photons, the mean photon number is expected to be 
(n) = \p2 and n' = 3^/2 in this case. 

Again, we carried out simulations which confirmed the 
existence of the expected resonance. This time, however, 
the simulation record looks very different from the results 
of the previous section, qs shown on Fig. [4j The mean 
photon number is low, and so is the average transmitted 
flux, but every time a photon is detected, we now see a 
drastic increase in the mean photon number in the cav- 
ity. This has a simple explanation, since the state prior 
to the detection event is a superposition of the ground 
state of the system and a doubly excited dressed state. 



The back action of the detection of a single photon, the 
quantum jump, is implemented by the action of the field 
annihilation operator a, which replaces the state before 
the detection event by 

|*U oc a {|G, 0) + v^(|G, 2) ± \E, 1))} 

oc y/2\G,l) ± \E,Q). (6) 

The mean photon number in this state is 2/3, which is 
much larger than the potentially nearly infinitesimal time 
averaged photon number in the cavity. The figure shows 
the mean photon number evaluated with the stochastic 
wave function method, and the transient peaks confirm 
our simple analysis. The detection record also shows a 
strong bunching effect: a majority fraction of the jumps 
are followed by a second jump within the cavity lifetime. 
A dedicated experiment should be able to verify these 
non-classical intensity correlations, and in a future per- 
spective, one may even imagine the possibility to apply 
feed-back and modify the driving field immediately after 
the first click event and thus prepare a variety of other 
quantum states of the system [39] . 

C. Three-photon resonance 

More generally, the spectrum Eq. [4] predicts the exis- 
tence of n-photon resonances at the driving frequencies 
uj n = (nu c ± g^\fn)jn = uo c ± # e ff/ f y/n. For the specific 
experimental setup considered here, these resonances are 
well distinguishable, even for n > 2 : for example, the 
3-photon resonance is well separated by a few line widths 
from the 2- and 4-photon resonances. Of course, the 3- 
photon excitation is a higher order process with interme- 
diate non-resonant 1- and 2-photon virtual excitations of 
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FIG. 2: Dressed ladder states of the atoms-cavity system: a) without inter-atomic interaction, the system would be formally 
equivalent to two coupled oscillators and transition frequencies between the multiplicities are degenerate; b) in case of interaction 
between excited atoms strong enough for the atomic sample to be fully dipole-blocked, the system is described by the JCM, 
and the nonlinearity prevents simultaneous excitation of several multiplicities; c) for larger atomic samples or weaker dipole 
interaction, several Rydberg "bubbles" can appear, and the system's spectrum is described by " combinatorially dressed states", 
still preserving a nonlinearity for small enough bubble numbers. 
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FIG. 3: Detection record for one-photon resonance, a — 2tt x 
0.15 MHz. 



FIG. 4: Detection record for two-photon resonance, a 
2tt x 1.5 MHz. 



the system, and it competes with the non-resonant exci- 
tation of the system by lower order processes. All these 
processes are included in the full numerical simulations 
of the system. 

Fig. [5] shows a simulated detection record for a driv- 
ing field resonant with the triply excited dressed state 
of the atoms+cavity system. This time, our simple 
analysis suggests that clicks cause a quantum jump of 
the state from a superposition of the ground state and 



(|G,3 > ±\E,2 >) towards y/3\G,2 > ±y/2\E,l >, 
which has an average photon number of 8/5. Within 
the cavity lifetime, one should thus expect bursts of 2-3 
photon detection events. 

Such bunching of detection events is clearly seen in Fig. 
5 But when we zoom in on the detection record (see Fig. 
6), we observe somewhat surprising "bursts" of 3, 4 and 
even up to 7 detection events, which contradict our sim- 
ple picture of the dynamics. Looking more closely into 
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FIG. 5: Detection record for three-photon resonance, a 
2tt x 1.25 MHz. 




FIG. 6: Detection record for three-photon resonance (magni- 
fied view). Each click is indicated by a bullet symbol on the 
time axis, a — 2tt x 2 MHz. 



IV. COMBINATORIAL EXTENSION OF 
DIPOLE BLOCKADE MECHANISM 

In the previous sections, we suggested to take advan- 
tage of the full Rydberg blockade in a small atomic en- 
semble to effectively implement the JCM in the strong 
coupling regime and we studied the transmission proper- 
ties of such an ensemble-cavity system. In this section, 
we shall see how the coupling of the atoms to the cavity 
can actually be used to extend the range of the blockade 
mechanism beyond the limit fixed by the typical range 
£r ~ few jam of Rydberg dipole-dipole interaction. 

In the absence of interactomic interactions, i.e. for en- 
sembles with interatomic separations larger than £r, the 
atomic ensemble and the cavity mode behave like two 
coupled oscillators leading to two atom-field eigenmodes 
with frequencies uj c db g e ff (Fig. [5] a). Driving the cavity 
with a classical probe field tuned on one of these fre- 
quencies, one resonantly excites a coherent state of the 
corresponding eigenmode equivalent to product of coher- 
ent states of the field and the collective atomic oscillator 
comprising all the Rydberg excitation number states. In 
the opposite limit, i.e. for ensembles whose linear dimen- 
sion is of the order of £r, the dipole-dipole interactions 
between Rydberg atoms are so strong that they shift the 
multiply Rydberg excited states out of resonance: the en- 
semble is equivalent to a two-level atom, and the entire 
system is described by the JCM, Eq.^ (Fig. |2]b). 

Let us now turn to the intermediate regime where the 
extent of the ensemble is larger than £r but the Ryd- 
berg dipole-dipole interactions still does play an impor- 
tant role in the system. In that case, the interaction does 
not shift all the multiply excited states out of resonance 
and we may identify an unshifted doubly excited collec- 
tive state 



the state conditioned on the subsequent detection events, 
we have identified two mechanisms which seem to con- 
tribute to the observed effect: i) the system does not only 
populate the ground and triply excited dressed states, a 
small off-resonant excitation of higher lying states also 
occurs. The corresponding population is then amplified 
by the action of the annihilation operators and the y/n 
factors associated with the detections of the first, second 
and third photon, ii) the system evolves between clicks, 
both due to the driving field and due to the states of 
the form \/3|G, 2 > ±y/2\E, 1 > resulting from the first 
click not being eigenstates of the atoms-cavity coupling. 
Also, the small probability amplitude in excited states, 
which is being magnified by a quantum jump, contributes 
coherently to the build-up of further excited state popu- 
lation by the coherent probe field. We have qualitatively 
verified the last point by setting the probe amplitude a to 
zero after the first click in our simulations and observed 
a reduction of the number of transient detection events. 

This multi-photon burst phenomenon constitutes an 
interesting experimental effect, but we will not investi- 
gate it further at this point. 



\E 2 ) 



A 



\ e j, e k) {9j,9k\ \G) (7) 



K \rj-r k \>£ R 



in which the excited atoms (j, k) are always too far apart 
from each other to interact strongly. The normalization 
constant A can be evaluated by counting the number of 
states \ej,ek) which contribute to \E 2 ), i-e. the number 
of pairs of atoms in the ensemble which are separated 
by a distance larger than £r. This state is thus reso- 
nantly accessible from the first excited state \E) via the 
two-photon absorption from the cavity and the blue laser 
classical field. The Hamiltonian Eq. ([3| must therefore 
be complemented by the term 

2faj c \E 2 ) (E 2 \ + hg ef f (E 2 S+ e) \E 2 ) (E\ a + h.c. 



More generally, due to the size of the sample, the res- 
onant coupling by the cavity mode and laser fields to 
higher multiply excited states \Es) , \E±) , . . . in which 
n e = 3,4,... Rydberg atoms are separated by more than 
£r, may be possible. We can visualize the system as if 
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the Rydberg blockade phenomenon decomposes the en- 
semble into a number of "bubbles" of radius £r [40]. In 
each bubble only a single Rydberg excitation is possible 
and the number of atoms which can be simultaneously 
excited in the ensemble coincides with the number of 
bubbles, ~ -j^w- For a mesoscopic sample which can 
thus accommodate n& Rydberg excitations, the Hamilto- 
nian Eq. ([3| must therefore be complemented by [ni — 1) 
terms coupling \E ) = \G) to \E X ) = \E), |£a) to \E 2 ), 
\E 2 ) to \E S ), ... , |S nb _i) to \E nh ). 

The ensemble-cavity dynamics is thus governed by the 
Hamiltonian 



H = hLJ c a ] a + huj c \ Y^k\E k ) {E k \j 

+ a(^g k \E k )(E k _ 1 \^ + h.c. (8) 

The coupling strengths g k can be easily evaluated in 
the "bubble" picture, since each bubble provides an effec- 
tive two- level atom and the states \E k ) are the symmet- 
rically excited states of these effective two-level systems 
with precisely k excitations. These states, in turn, are 
equivalent to effective angular momentum states |j, m) 
with the quantum numbers j = and m = — 77,5/2 + k. 

The coupling strengths g k are given by the well-known 
matrix elements of the angular momentum raising opera- 
tor multiplied by the effective two-state coupling of each 
bubble, (fcl^lfc — 1), and we get 



9k 




(9) 



Setting k = 1, we recover the expected collective cou- 
pling g e fj to the first excited states, while 
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n b A 



(10) 



which vanishes for = 1 as it should, and which ap- 
proaches the expected ^/2g e ff in the non-blocked oscil- 
lator limit of large n^. 

Any finite value, and in practice a not too large value, 
of 77,5 thus breaks the coupled oscillator picture, and hence 
the excitation spectrum of the coupled system will not 
constitute an equidistant ladder of states. In particular, 
the resonant excitation of the first dressed states with 
a single Rydberg atom component is not resonant with 
the next excitation step, and the system thus behaves 
as if the Rydberg blockade extends over the entire sys- 
tem across different "bubbles" (Fig. [5]c). We recall that 
the absence of simultaneous excitation of several bubbles 
is not due to their mutual interaction but due to the 
coupling strength to higher excited states which is mod- 
ified because of the number of pairs that contribute to 



this excitation. The blockade behavior is thus due to a 
combinatorial effect. Diagonalizing the Hamiltonian (J8|, 
within the space of two excitations, we find the eigenval- 
ues 0, ±w2 — -^h\/2g e ff, and hence the field resonant 



with the first dressed state is detuned from the 

two-photon resonance by 2 ^1 — ^Jl 



2n b ^ e ff ) 



2n b 



This analysis suggests that we can only extend the 
blockade over a few, say less than ten, "bubbles". We 
intend to perform more detailed numerical studies of the 
collectively excited atoms. The picture with a given fixed 
set of bubbles of atoms is not an exact one as any atom 
has the same amplitude to be excited and hence form 
the center of a bubble of excluded excitations around it, 
and the full many body state is a superposition of cor- 
respondingly different bubble configurations [41-43 . We 
do expect, however, that our combinatorial argument is 
robust to finer details of the description and that any 
two-atom collective component will not be excited at the 
same frequency as the dressed states with a single collec- 
tive excitation. 



V. DISCUSSION 

We have suggested to use the collective coupling of 
~ 1000 atoms to a cavity field combined with the Ryd- 
berg blockade effect to restrict the atomic ensemble to an 
effective two-level system. Physical parameters realisti- 
cally reach the JCM of strong coupling CQED, and we 
have shown that exploration of higher excited states of 
the dressed states ladder of this model should be possible 
by optical transmission experiments. 

In a longer perspective, many more possibilities can be 
investigated. One may take advantage of the atomic sub- 
level structure to significantly enrich the phenomenology 
with collective qubit encoding [44], multi-mode atomic 
storage states, STIRAP processes |45], and several quan- 
tum field components may couple to Rydberg states with 
different pairwise interaction properties [46 to accommo- 
date effective optical non-linearities [47] . 

In this work, we disregarded the role of atomic mo- 
tion, but we note that even though the blockade regime 
precludes strong interatomic forces as only single atoms 
are excited, the position-dependent collective field-atom 
coupling may induce mechanical motion entangled with 
the field and atomic internal state degrees of freedom. If 
the atoms form a Bose-Einstein condensate of interacting 
atoms, as in [17], collective excitations will then act as 
micro- mechanical degrees of freedom [19] [20] , with po- 
tentially strong coupling to the collective internal state 
two-level system. 

Another interesting direction of research, involves the 
examination of larger systems. In such systems, the Ryd- 
berg blockade is only efficient in blockade spheres around 
each atom, and remote atoms may be simultaneously ex- 
cited. We have proposed that the combinatorics of the 
atoms-cavity coupling may distort the ladder of energies 
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of singly and multiply occupied Rydberg states, and thus 
maintain the blockade effect over a larger atomic ensem- 
ble. This should be more quantitatively addressed in 
future work. 



mon Michels. 
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